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Orbital eccentricity of the symbiotic star MWC 560 * 
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We present projected rotational velocity measurements of the red giant in the symbiotic star MWC 560, using the high- 
resolution spectroscopic observations with the FEROS spectrograph. We find that the projected rotational velocity of the 
red giant is v sin i— 8.2 ±1.5 km s , and estimate its rotational period to be P TO t = 144 - 306 days. Using the theoretical 
predictions of tidal interaction and pseudosynchronization, we estimate the orbital eccentricity e = 0.68 — 0.82. We briefly 
discuss the connection of our results with the photometric variability of the object. 
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1 Introduction 

MWC 560 (V694 Mon) was discovered as an object with 
bright hydrogen lines (Merrill & Burwell 1943). It is a sym- 
biotic binary system, which consists of a red giant and a 
white dwarf (Tomov et al. 1990; Michalitsianos et al. 1993). 
The orbital period is supposed to be P or b ~ 5.3 yr (Doro- 
shenko, Goranskij & Efimov 1993). The most spectacular 
features of this object are the collimated ejections of matter 
with velocities of up to ~ 6000 kms -1 (Tomov et al. 1992; 
Stute & Sahai 2009) and the resemblance of its emission 
line spectrum to that of the low-redshift quasars (Zamanov 
& Marziani 2002). 

The jet ejections are along the line of sight and the sys- 
tem is seen almost pole-on (i < 16°). This makes it diffi- 
cult to obtain the orbital parameters using the radial velocity 
variations, and impossible to observe effects such as eclipse 
or illumination. 

To improve our understanding of this object, we aim 
here to: (1) measure the projected rotational velocity of the 
red giant (v sin i) from high resolution spectra, and (2) use 
our measured v sin ito find the eccentricity e of the orbit, on 
the basis of the theory for pseudosynchronization in binary 
stars. 

2 Spectral Data and v sin i measurement 
2.1 Observations 

We analyzed 21 high resolution spectra of MWC 560, ob- 
tained with the FEROS spectrograph. FEROS is a fiber fed 
echelle spectrograph, providing high resolution of A/ A A = 

* based on: observations obtained in ESO program 074.D-01 14 
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48000 and an wide wavelength coverage from about 4000 A 
to 8900 A in one exposure (Kauferet al. 1999). 16 of the an- 
alyzed spectra are from our observations of the MWC 560 
with the 2.2-m telescope at ESO, La Silla (ESO program 
074.D-0114). The remaining spectra (5), are downloaded 
from FEROS Spectroscopic Database at the LSW Heidel- 
berg and were obtained when FEROS was mounted on the 
1 .52-m telescope at La Silla. A log of observations is given 
in Table Q] Details of the data processing are given in Za- 
manov et al. (2007). In Fig.Q]we show the average spectrum 
derived from the 16 spectra obtained by 2.2-m telescope. 

Symbiotic stars have composite spectra with three main 
sources of radiation - red giant, hot component and nebula. 
Modeling of a few objects is given in Skopal (2005). Fig.Q] 
also shows a synthetic spectrum of a M5III star with param- 
eters: T cff = 3424 K, log g = 0.5, £ = 3 kms -1 , v s'mi= 
8.2 km s~\ and 64% contribution of the giant at A = 8800 A 
(see 2.2.). 

2.2 v sin i measurement 

We measure v smi by comparing the FWHM of spectral 
lines of observed and synthetic spectra. The procedure is 
similar to that described in Fekel (1997). 

We synthesized spectra by using the code SYNSPEC 
(Hubeny, Lanz & Jeffery 1994) in the spectral region A = 
8750 - 8850 A. Atmospheric parameters T cff = 3424 K, 
log<? = 0.5 (typical for a M5III star, see van Belle et al. 
1999), the instrumental profile for our FEROS spectra, and 
solar abundances were the input parameters for the code. 
LTE model atmospheres were extracted from Kurucz's grid 
(1993). The VALD atomic line database (Kupka et al. 1999) 
was used to create a line list for spectrum synthesis. 



www.lsw.uni-heidelberg.de/projects/instrumentation/Feros/ferosDB/ 
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Table 1 Spectral observations of MWC 560. The columns 
in the table represent: (1) modified Julian Date (MJD) of 
the start of the exposure, (2) exposure time, and the (3) 
measured projected rotational velocity (v sin i) of the mass 
donor. 



MJD 


exp-time 


v ami 




[sec] 


[kms -1 ] 


2.2-m telescope 






53314.279 


300 


8.5 


53314.283 


300 


8.6 


53314.287 


300 


7.9 


53314.291 


300 


7.8 


nil a r\f\c 

53314.295 


300 


8.0 


53314.299 


300 


8.2 


53314.303 


300 


8.2 


53314.308 


300 


8.5 


53314.312 


300 


8.8 


53314.316 


300 


7.8 


53314.320 


300 


8.3 


53314.324 


300 


7.7 


53314.328 


300 


8.6 


53314.332 


300 


7.9 


53314.336 


300 


8.1 


53314.340 


300 


8.5 


1.52-m telescope 






51133.223 


1200 


7.8 


51131.248 


2400 


8.9 


51135.283 


1200 


7.8 


51136.259 


900 


7.9 


51141.315 


1200 


8.2 



We adopt a micro turbulent velocity £ = 3 kms -1 . This 
value is the same as accepted by Fekel, Hinkle & Joyce 
(2004) and it is close to the values for other symbiotics. 
Schmidt et al. (2006) calculated £ = 2.2 km s" 1 for CH Cyg.| 
Wallerstein et al.(2008) used £ = 2.2 km s" 1 for T CrB. The 
relation between £ and log g (Gratton, Carretta, & Castelli 
1996) gives £ k 2.1 kms -1 for our case. 

A grid of synthetic spectra for projected rotational ve- 
locities from kms -1 to 20 kms -1 was calculated, v sini 
was calculated by measuring the full width at half max- 
imum (FWHM) of twelve observed spectral lines. These 
were compared to the spectral line half- widths from the syn- 
thetic spectra. We thereby derive a mean v sin i= 8.2 km s _1 ,| 
and a standard deviation of the mean a — 0.4 km s -1 using 
all 21 spectra, we note that if we adopt a microturbulent ve- 
locity £ = 2 kms -1 we would have measured v sin«= 9.4 
km s _1 . For the error of our measurements we adopt a con- 
servative value ±1.5 kms , which also includes the uncer- 
tainty of £ = 3 ± 0.5 km s _1 and the error of our method. 
We therefore estimate our final measurement to be v sini= 
8.2 ±1.5 kms" 1 . 



3 Tidal interaction in MWC 560 

3.1 Parameters of MWC 560 

Meier et al. ( 1 996) classified the cool component of MWC 5 60| 
as M5IH-M6III giant. Miirset & Schmid ( 1 999) give M5 .5III| 
- M6III as its stellar type. According to van Belle et al. 
(1999), a typical M5III giant has the following parameters: 
R e = 139.6±5.6R ,T eff = 3424 K, L g - 2410 L Q , and a 
typical M6III giant: R s = 147.9 ± 7.7 R Q , T cS = 3375 K, 
L e ~ 2560 L G . We adopt for the M giant in MWC 560 
R g = 140 ± 7 R G and L g - 2400 L s . 

For the masses of the red giant and white dwarf we will 
assume M g = 1.7 M Q and M wc i = 0.65 Mq, respectively. 
These values are the average masses of the components of 
the symbiotic stars (Mikolajewska 2003). Concerning the 
orbital period, we use P or b = 1931±162 day (Gromadzki et 
al. 2007). With the above values of the parameters assumed, 
we derive the semimajor axis of the orbit to be a « 860 Rq. 

The jet (orbit) inclination angle to the line of sight is 
i < 16°, and the white dwarf accretes at a rate Af acc w 
5 x 10~ 7 M Q (Schmid et al. 2001). 

3.2 Synchronization and circularization time scales 

The physics of tidal synchronization for stars with convec- 
tive envelopes has been analyzed several times. There are 
some differences in the analysis of different authors, leading 
to varying synchronization timescales. We use the estimate 
from Zahn (1977, 1989). The synchronization timescale in 
terms of the period is 

'syn 



800 



Rl 



-Cb yr, (i) 



where M g and M w( j are the masses of the giant and white 
dwarf respectively in Solar units, and R g , L g are the radius 
and luminosity of the giant, also in Solar units. The orbital 
period P or b is measured in days. 

Following Hurley, Tout & Pols (2002), the circulariza- 
tion time scale is: 

/ i - \ 8 

92(1 + 92) 



1 



R„ 



21 (k 

^"circ 2 y T J \^ Q> 

where is the mass ratio — M W( x/M g 
is derived from Rasio et al. (1996): 

k\ 2 / C onv M311V _i 

j) ~ 21 w M g yI ' 
where R cnv is the depth of the convective envelope, M, 
is the envelope's mass, and 

M env R env (Rg 2-^env) 



(2) 

InEq.El (k/T) 
(3) 



0.4311 



3L e 



yr (4) 



is the eddy turnover time scale (the time scale on which 
the largest convective cells turnover). The numerical factor 

/conv IS 

\-[—^-) , (3) 



/ conv = min 



2r„ 
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Fig. 1 The average spectrum of MWC 560 in the range A = 8790 — 8850 A, obtained from 16 exposures with FEROS 
spectrograph at the 2.2-m telescope at La Silla, and a synthetic spectrum using T ff=3424 K, log g = 0.5 £ = 3 kms -1 , 
v sin i= 8.2 km s -1 , and 64% contribution from the giant. 



where P t id is the tidal pumping time scale given by 



1 



1 



-Porh 



l 

-Prot 



(6) 



For MWC 560 we assume for the red giant to have R cm = 
0.9 i? g and M env = 1.0 M (Herwig 2005). Using these 
parameters, we calculate P t id = 248 days, / conv = 1, 
r conv = 0.476 yr, and (k/T) = 0.12 vr" 1 . Then from 
Eq.[T]and Eq.|2] follows that the synchronization and circu- 
larization time scales are r sync = 2.6 x 10 4 yr and r c i rc = 
3.1 x 10 6 yr respectively. 

Following Hut (1981), we estimate the pseudosynchro- 
nization timescale r DS : 



7 



3(a - 3) 



(7) 



where a is a dimensionless quantity, representing the 
ratio of the orbital and rotational angular momentum: 



92 1 



a 

1 + 92 r 2 \R t 



(8) 



where r g is the gyration radius of the giant. For a red giant 
we adopt: r g w 0.3 (Claret, 2004, 2007). We calculate a = 
800 andr ps =9.1 x 10 3 yr. 



3.3 Lifetime of the symbiotic phase 

The typical lifetime of a symbiotic star is r ss ~ 10 5 yr (Yun- 
gelson, et al. 1995; Lii, Yungelson, & Han 2006). In the case 
of MWC 560, we can estimate from the rate of accretion 
on the white dwarf, Mace ~ 5 x 10~ 7 M Q (Schmid et al. 
2001), that it will take 10 6 yr to accrete ~ 0.5 M Q from the 
envelope of the red giant companion. Because the giant also 
losses mass via stellar wind, we find that the lifetime of the 
symbiotic phase of MWC 560 should be r ss ;$ 10 6 yr. 



For MWC 560 we have therefore the situation in which 
T ps < T syn < t ss < T c i rc . This means that the symbiotic 
phase is long enough that the tidal forces can synchronize 
(pseudosynchronize) the rotation of the red giant. On the 
other hand, the value of T c j rc demonstrates that the symbi- 
otic lifetime of MWC 560 is shorter than the circulariza- 
tion time, and therefore the orbit can be eccentric. This is in 
agreement with the observational evidences found by Fekel 
et al. (2007) that the symbiotic stars with P or b > 800 days 
tend to have eccentric orbits. 

Schmutz et al. (1994) estimated that the timescale for 
circularization in SY Mus is ~ 10 times longer then the 
synchronization time. For MWC 560, we find even higher 
ratio: r circ /r syn w 45. 

The above implies that in MWC 560, the red giant is 
probably more or less synchronized, but the orbit is not cir- 
cularized. 

3.4 Rotation and pseudosynchronization 

In a binary with a circular orbit the rotational period of the 
primary, P rot , reaches an equilibrium value at the orbital 
period, P rot = -Porb- However, in a binary with an eccen- 
tric orbit, the tidal force acts to synchronize the rotation of 
the mass donor with the motion of the compact object at 
the periastron - the effect called pseudosynchronous rota- 
tion (Hall 1986). The corresponding equilibrium (i.e. pseu- 
dosynchronization) is reached at a value of P rot which is 
less than P orD , the amount less being a function of the or- 
bital eccentricity e. Hut (1981) showed that the period of 
pseudosynchronization, P ps , is : 



Pv 



fl + 3e 2 



l )(l 



1 



45 „4 
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JL P 6 
16 c 



(9) 



When the eccentricity tends to zero, P ps tends to P or D . 
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3.5 Orbital eccentricity of MWC 560 

In order to determine the orbital eccentricity of MWC 560, 
we first need to calculate P rot for the mass donor. We use 

P ro t = (10) 
v sin i 

The underlying assumption is that the rotational axis of the 
mass donor is perpendicular to the orbital plane. 

Using the value for v sini= 8.2 ± 1.5 kms" 1 (SectfS]), 
i? g = 140 ± 7 R Q , and i = 12° - 16° (Schmid et al. 2001), 
we calculate P rot = 144 — 306 days. This value is consid- 
erably less than the orbital period, P or b = 1931 ± 162 days. 
Following the results in Sect l3.21 MWC 560 should be close 
to synchronization or pseudosynchronization, and P rot = 
P ps . Finally, using Eq. [9] we can therefore estimate the or- 
bital eccentricity to be e = 0.68 — 0.82. 

4 Discussion 

In most of the symbiotics the activity of the hot components 
is irregular/aperiodic. However, in MWC 560 it has a pe- 
riodic character. An analysis of the historical light curve 
(Luthardt 1991) and B-band photometry, as well as the pe- 
riodogram analysis of V band and near-IR observations re- 
veal that this periodicity has remained in-phase for over a 
century (Doroshenko et al. 1993; Gromadzki et al. 2007). 
The most natural explanation for this phenomenon is or- 
bital modulation: around periastron, the mass accretion rate 
increases, which causes an increase of the accretion lumi- 
nosity and the brightness of the accretion disk, which pro- 
duces the observed modulation of the light curves (see also 
Doroshenko et al. 1993; Gromadzki et al. 2007). 

With the parameters assumed above, we can calculate 
the distance between the components at the periastron: r p = 
a(l - e) = 180 - 230 R Q . It follows that the red giant fills 
the Roche lobe for about 5% - 25 % of the orbital period. 
This is in agreement with the suppositions of Gromadzki 
et al. (2007), that the orbital eccentricity and corresponding 
Roche lobe overflow at the periastron is the reason for the 
observed photometric variability. 

It has been noted that the rotational period of the mass 
donor is considerably shorter than the orbital period (Za- 
manov et al. 2008) in a few jet-ejecting symbiotics. Our 
findings here pose the question whether their orbits are also 
eccentric. 



Conclusions 

In this note, we presented new measurements of the pro- 
jected rotational velocity of the red giant in the symbiotic 
star MWC 560. We find that v sini= 8.2 ± 1.5 kms" 1 . 
On the basis of the theory of tidal interaction in binaries, 
we calculate that the orbit should be highly eccentric, with 
e £ 0.7. Our findings support the model that the observed 
photometric variability of MWC 560 is connected with high 
orbital eccentricity and Roche lobe overflow at periastron. 
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